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In this paper, we have investigated theoretically the influence of atomic collisions on the be- 
haviour of a one-dimensional Bose-Einstein condensate inside a driven optical cavity. We develop the 
discrete-mode approximation for the condensate taking into account the interband transitions due 
to the s-wave scattering interaction. We show that in the Bogoliubov approximation the atom-atom 
interaction shifts the energies of the excited modes and also plays the role of an optical parametric 
amplifier for the Bogoliubov side mode which can affect its normal-mode splitting behaviour. On 
the other hand due to the atomic collisions the resonance frequency of the cavity is shifted which 
leads to the decrease of the number of cavity photons and the depletion of the Bogoliubov mode. 
Besides, it reduces the effective atom-photon coupling parameter which consequently leads to the 
decrease of the entanglement between the Bogoliubov mode and the optical field. 

PACS numbers: 03.75.Gg, 42.50.Wk, 67.85. Hj, 37.30. +i, 03.67.Bg 



I. INTRODUCTION 

The emerging field of optomechanics concerns with the 
study of the mechanical effects of light on mesoscopic and 
macroscopic mechanical oscillators. This phenomenon 
has been realized in the optomechanical systems con- 
sisting of an optical cavity with a movable end-mirror 
or with a membrane-in-the middle. The radiation pres- 
sure exerted by the light inside the optical cavity cou- 
ples the moving mirror or the membrane which acts as a 
mechanical oscillator to the optical field. This optome- 
chanical coupling has been employed for a wide range of 
applications such as the cavity cooling of microlevers and 
nanomechanical resonators to their quantum mechanical 
ground state PUS] , producing high precision detectors for 
measuring weak forces and small displacements and also 
providing a good approach for fundamental studies of the 
transition between the quantum and the classical world 

[SHE]. 

On the other hand an alternative path to the studies 
of cavity optomechanics has been provided experimen- 
tally by systems consisting of ultracold atomic ensembles 
trapped in optical cavities [T0HT3] where the excitation 
of a collective mode of the cold gas plays the role of the 
vibrational mode of the mirror. The standing electro- 
magnetic wave inside the cavity forms a periodic poten- 
tial, the so-called optical lattice, in which the cold atoms 
exhibit phenomena typical of solid state physics like the 
formation of energy bands and Bloch oscillations |14) . 

In such systems with high finesse cavities the atom- 
light interaction is enhanced because the atoms are col- 
lectively coupled to the same optical mode. Besides, in 
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the dispersive regime where the laser pump is far detuned 
from the atomic resonance the excited electronic state 
of the atoms can be adiabatically eliminated and con- 
sequently the only degrees of freedom of atoms will be 
their mechanical motions [TSHTT] . For low photon num- 
bers when the optical grating produced by the intracav- 
ity optical field is very shallow, one can approximately 
restrict the dynamics to the first two motional modes 
. In more recent theoretical investigations [23] US] 
it has been shown that the simple two-mode model of 
a Bose-Einstein condensate (BEC) can be improved by 
considering higher motional modes. 

In spite of similarities between the two kinds of op- 
tomechanical systems (with a moving mirror and with a 
BEC), there are some essential differences between them. 
Firstly, in contrast to the moving mirror of the optome- 
chanical sytems which is placed in the harmonic poten- 
tial of its spring, the excitation modes of the BEC are 
not based on the presence of such an external harmonic 
potential |12) . Secondly, their parameters are realized in 
different regimes; the frequency of the excitation mode 
of the BEC, i.e., recoil frequency, is well below the dipole 
coupling strength while the oscillation frequency of a 
moving mirror is of the same order of magnitude as the 
coupling strength (T3] [25] . On the other hand in addi- 
tion to pumping the cavity from one of the end mirrors, 
in the cavities equipped with atomic gas it is also possible 
to pump the atoms from the side of the cavity [27]. 

One of the most important characteristics of the many- 
body systems is the two-body collision which can affect 
the properties of the system. So in order to study the 
dynamics of a BEC gas in a realistic experimental sit- 
uation, it is necessary to take it into account. In the 
optomechanical systems containing the atomic gas there 
are two kinds of atom-atom interactions which are the 
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origins of the system nonlinearities. Firstly, due to the 
atom-photon interaction the potential acting on the con- 
densate depends in a highly nonlocal and nonlinear way 
on the condensate itself which leads to the long range 
atomic interaction mediated by the cavity field [23 OH] . 
On the other hand there is an intrinsic nonlinearity due 
to the s-wave scattering which can take place at arbi- 
trary momentum values and causes a broadening of the 
atomic momentum distribution [2D] due to the intraband 
transitions and also can scatter atoms to the other bands 
(interband transitions) O HDJ • Furthermore, there is an- 
other kind of nonlincaritiy due to the finiteness of the 
particle number of the system which can be neglected 
in the thermodynamic limit where the total number of 
atoms is very large [ID] . 

In this work, we are going to extend the two-mode 
model considering the effects of atom-atom interaction 
for a one-dimensional BEC in an optical cavity. We con- 
sider the system in the low photon regime but will not 
restrict our treatment to the weak atom-atom interac- 
tion. In fact the atomic collisions can populate several 
nonzero quasimomenta of the energy bands (intraband 
transition) and can also cause the atoms to be scattered 
to the other bands (interband transitions). In order to 
have a simplified optomechanical model we take into ac- 
count just the lowest nozero quasimomenta in the first 
Brillouin zone . We show that in the Bogoliubov approx- 
imation the atom-atom interaction not only shifts the 
energies of the excited modes but also play the role of 
an oprtical parametric amplifier (OPA) for the Bogoli- 
ubov side mode which can affect its normal-mode split- 
ting behaviour [DJ . On the other hand due to the atomic 
collisions the resonance frequency of the cavity is shifted 
which leads to the decrease of the number of cavity pho- 
tons and the depletion of the Bogoliubov mode. These 
results are in good agreement with those obtained by 
the numerical solutions of the full description of Gross- 
Pitaevskii Equation (GPE) [551 \M\. 

The paper is structured as follows. In Sec. II we will 
first give a thorough theoretical description of the many- 
body system under consideration and then derive a sim- 
plified optomechanical model. In Sec. Ill the quantum 
Langevin equations (QLEs) are derived and linearized 
around the semiclassical steady state. In Sec. IV we 
study the mean-field solutions and fluctuations of the 
system and in Sec. V we derive the spectrum of the 
Bogoliubov mode and investigate the normal mode split- 
ting (NMS). Finally, our conclusions are summarised in 
Sec. VI. 



II. THEORETICAL DESCRIPTION OF THE 
SYSTEM 

We are going to study a gas of N ultracold bosonic 
two-level atoms with mass M and transition frequency 
Lu a inside the optical lattice of a single-mode, high-finesse 
Fabry-Perot cavity with length L. The cavity is driven at 
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FIG. 1. (Color online) N two-level atoms trapped in an opti- 
cal cavity interacting dispersively with a single cavity mode. 
The cavity mode is driven by a laser at rate r\ and the decay 
rate is k. 



rate r\ through one of its mirrors by a laser with frequency 
LJp, and wavenumber K = 2n/X = lj p /c. We assume 
the BEC is confined in a cylindrically symmetric trap 
with a transverse trapping frequency uj± and negligible 
longitudinal confinement along the x direction (Fig[I]). 
In this way we can describe the dynamics within an ef- 
fective one-dimensional model by quantizing the atomic 
motional degree of freedom along the x axis only. 



A. The General Form of The Hamiltonian of a 
BEC Inside an Optical Lattice 

In the dispersive regime where the laser pump is far de- 
tuned from the atomic resonance (A„ = uj p — uj a exceeds 
the atomic linewidth 7 by orders of magnitude), the ex- 
cited electronic state of the atoms can be adiabatically 
eliminated and spontaneous emission can be neglected 
|15j . In the frame rotating at the pump frequency, the 
many-body Hamiltonian reads 

[L/2 

H = -ftA c a f a+itvq{a) -a)+ / dx& (x)H Q ^(x)+H aa , 

J -L/2 

(1) 

where a is the annihilation operator for a cavity photon 
and A c = lo p — lo c is the cavity-pump detuning. Hq is the 
single-particle Hamiltonian of an atom inside the optical 
lattice of the cavity and H aa is the atom-atom interaction 
that are respectively given by 

H = p 2 /2M + hU a cos 2 (Kx)a?a, (2a) 
1 f L/2 

H aa = -U s / dxtf(x)&(x)9(x)*(x). (2b) 

1 J-L/2 
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Here Uo = g 2 /A a is the optical lattice barrier height 
per photon and represents the atomic backaction on the 
field, go is the vacuum Rabi frequency, U s — 47r (L a3 and 
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a s is the two-body s-wave scattering length [TBI 02] ■ The 
second term in the Hamiltonian of Eq. (|2^) is a periodic 
potential of period A/2. 

It is well-known that the eigenf unctions of a particle 
inside a periodic potential are the Bloch functions tpu^ix) 
with eigenvalues e„ jg where v is the Bloch band index and 
q G [-2ir/\, 2ir/\] is the quasimomentum of the particle 
|18j . If there are / periods inside the cavity, then L = 
IX/2. Using the Born- Von Karman periodic boundary 
condition %p v . q (x + L) = ip Uiq (x) and based on the Bloch 
theorem, 



tp u , q (x + IX/2) 



(3) 



it is deduced that qlX/2 = 2rrnr, where (m G Z). In 
this way the quasimomentum q m in the first Brillouin 
zone can be written in terms of optical wave number, 
i.e., q m — 2mK/l where —1/2 sC m ^ 1/2. Now using 
the Bloch theorem, the eigenfunctions ipu,q m (x) can be 
written as 



1 



°u u (x). 



(4) 



Here u v {x) is a periodic function with the optical lattice 
period (A/2) and reads 



u v {x) 




i2nKx 



(5) 



By substituting Eq.Q into Eq.Q ,the eigenfunctions of 
the single-particle Hamiltonian Hq is obtained as 



i{q m +2nK)x 



(6) 



The second quantized atomic wave field ^(x) can be 
expanded as a Fourier series in terms of the Bloch func- 
tions, i.e., 



= ^2 K,q m i>u,q m {x). 



(7) 



Here b v ^ qm {b\, qm ) is the annihilation (creation) operator 
for the atomic field that annihilates (creates) a parti- 
cle in a state determined with the Bloch band index v 
and quasimomentum q m . By using this expansion the 
Hamitlonian of the system is diagonalized and one can 
investigate the Bloch band structure of the condensate 
[Jj5]. However, the Fourier coefficients c v>n can only be 
determined numerically and there is no closed form for 
the Bloch functions. Instead, by substituting Eq.Q into 
Eq. ([7]) one can obtain a Fourier expansion of the atomic 
field in therms of plane waves 



J2K(n+m/l)x 



(8) 



In the derivation of Eq.Q we have done a Bogoliubov 
transformation 



(9) 



where &„ j?n (w n m ) is the atomic field operator that an- 
nihilates (creates) a particle in a state determined with 
the band index n and quasimomentum q m . 

By substituting Eq.Q into Eq.Q the Hamiltonian of 
the systems can be written as 

H = -hA^a + ihr){a) - a) + e nm b\ m b nm 

+ - hU a f a b\ m (6 n _ 1 , m + b n +x,m + 2&„, m ) 

+H aa . (10) 
In these sums n G Z, —1/2 ^ m ^ 1/2, and 



{ n + j) 



(11) 



are the energy eigenvalues in which Er — h 2 K 2 /2M is 
the recoil energy. As seen, the field expansion of Eq.([8]) 
can no longer diagonalize the Hamiltonian of the system. 
In fact that part of Hamiltonian which corresponds to 
the atom-photon interaction remains nondiagonal. Due 
to the interaction with the optical filed, the atoms are 
scattered to the nearest bands (An = ±1) while their 
quasimomentum (m) remains unchanged. 



B. Optomechanical Model 

Since the optical potential is symmetric with respect to 
the origin, the Hamiltonian of the system has the parity 
symmetry. So if the sytem starts from a homogeneous 
BEC, only the cosine parts of the exponential functions 
in Eq. |8]) are excited because of the parity conservation. 
In the limit of weak photon-atom coupling, when either 
the photon number inside the cavity or Uq is small, the 
lowest bands n = ±1 can be excited by fluctuations re- 
sulting from the atom- light interaction [55] . On the other 
hand, the s-wave scattering populates fluctuations with 
arbitrary n and m |20j . In a very simplified optomechan- 
ical model one can consider the scattering from m = 
to m = ±1 in the lowest band n = 0. In this way the 
atomic field operator (Eq.®) can be truncated as 



*(x) = -^j- c oo + y ^cio cos(2if a;) + y ~c i cos(2Kx/l), 

where we have done the following Bogoliubov transfor- 
mations 

1 



V2 



(13a) 
(13b) 
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In the case that the system does not have parity sym- 
metry, for example when the BEC is inside a ring cavity, 
one should also consider sine modes with annihilation 
operators 



0, 



(14) 



which in our model have been set aside 2TJ [22 . By 
substituting the atomic field operator, Eq.(12|, into the 
Hamiltonian of Eq.Q, one can arrive at the following 
Hamiltonian 

H = -h/\ c a)a + ihrj(a - a f ) + 4E R (c\ c 10 + ^4i c oi) 



(15) 



where 



H, 



„, - -hU o) a(c 00 c o + c} cio + c 01 c i 



„t 



l t 1 

7f c °° Cl ° + 7f CooC i°' ) ' 



(16a) 



H aa — ^^j~^( c Oo( C 10 + C 0l) + c Oo( C 10 + C 0l) 



+4t<%i + c io c oi + Mo c io4i c oi 



J2 2 



J2 2 , +2 2 



+ 2 C oi c oi + 2 C io c io + c oo c ooy 



(16b) 



The first two terms in Eq.(15) denote the cavity and 



the pump Hamiltonians. The third term is the energy of 
the side-modes Cio and coi which is just the third term 
in Eq.( 10 1 with energy eigenvalues e nm given by Eq.( 11 1. 



The Hamiltonian H ac denotes the atom-photon interac- 
tion and H aa is the atom-atom interaction. Furthermore, 
uj sw — 8irha s N / M Lw 2 is the s-wave scattering frequency 
and w is the waist of the optical potential. 

This Hamiltonian can be further simplified since for 
weak optical fields and large N the depletion of the ini- 
tial condensate remains weak. So we can use the Bogoli- 
ubov approximation and treat the zero-momentum m ode 
classically, i.e., c 00 — > vN (211 [29]. In this way Eqs.(15) 
reduce to the following forms 

H = —hK c a) a + ihrj(a — a?) + h(4uj]i + oj sw )c 1q Cio 



Una = 



4 , 

-hU a^a(c w + c\ ) 



+ ^hU a a^ a(c\ Q c 10 + cJiCoi)) 



H aa — -7fu*} sw (C$ + c\ + Cqj + C J Q1 ) . 



J2\ 



(17a) 

(17b) 
(17c) 



Here A c = A c — NUq/2 is the effective Stark-shifted de- 
tuning. As seen from the Hamiltonian of Eq. (JTrjb) there 
are two kinds of optomechanical coupling. The first term 



of Eq.(17j) is the linear radiation pressure which cou- 



ples the side mode Cxo to the optical field with the op- 
tomechanical coupling constant ^j^-hUo while the sec- 
ond term is the nonlinear optomechanical coupling of 
the two side modes cio and Coi with the optical field. 
In the atom-atom interaction Hamiltonian H aa we have 
disregarded all terms proportional to u) sw /N. The influ- 
ence of s-wave scattering has partly appeared as a shift 
in t he si de modes energies (the third and fourth tems in 



Eq.(17a)). Since hu sw = 2U s tia, where ha = N/Lw is 



the density of atoms, it is apparent that the energy shits 
obtained here is the same as that of Bogoliubov theory 
|25j . On the other hand the role of H aa in Eq.(17l for 



the Bogoliubov side modes is very similar to that of an 
optical parametric amplifier (OPA) in an optomechani- 
cal system with the nonlinear gain parameter uj sw [9] . In 
Sec. V we will show how this parameter can affect the 
NMS behaviour of the coupled Bogoliubov mode and the 
cavity field. 



III. DYANAMICS OF THE 
OPTOMECHANICAL SYSTEM 



The dynamics of the optomechanical system described 
by Eqs.(17) is fully characterized by the following set of 
written in the frame rotating at the 



nonlinear QLEs 
input laser frequency, 



• ( ' a \ iU r V2N . t . t 
a = (iA c - n)a — a l—^— ( c io + cio) + c[ c w 



+4i c oi] -v + €, 



(18a) 



cio = -(Mdio +7) c io - i—j— U a'a —a<ac w 



~UJswC w + flO, 



a i \ iU o+ i t , t 

coi = -{iuqi +7)coi — a'ac i - -u sw c' n + ,/, 



where cjiq = 4u)r + u sw and wqi = jsUJr + u s 



(18b) 

01- 

(18c) 
Here 

k and 7 characterize the dissipation of the cavity field 
and collective density excitations of the BEC, respec- 
tively. The cavity-field quantum vacuum fluctuation 
£(t) satisfies the Markovian correlation functions, i.e., 
- (n ph + l)S(t-t'), (ZHt')Z(t)))=n P h5(t-t') 
with the average thermal photon number n p h which is 
nearly zero at optical frequencies [3T]. Besides, fw{t) 
and foi{t) are the thermal noise inputs for the two side 
modes of BEC which also satisfy the same Markovian 
correlation functions as those of the optical noise. The 
noise sources are assumed uncorrelated for the different 
modes of both the matter and light fields. 
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A. Linearization 

Now we are going to study the weak excitations of 
the condensate from its ground state. Such excitations 
include small deviations of both the atomic wave func- 
tion and the optical field from their respective stationary 
states. So we decompose each operator in Eqs.(18l as 



the sum of its steady-state value and a small fluctuation. 
By substituting a = a + Sa , cm = yN j3\ + 5c\q and 
cqi = VN 00 + Scqi into Eqs.( 18 ) one can obtain a set of 



nonlinear algebraic equations for the steady-state values, 



3 ' 



Pi = 



V2 U a 2 



^ )2 + 7 2 



0o = 0, 



(19a) 
(19b) 
(19c) 



where we have assumed a, (3q and /3i are real numbers 
[23] and A d = A c - \NU fi\{fi\ + \/2) is the effective 

detuning, Cl^' — lj w ± \lo sw and lo w = w 10 + \\Jqc? . 
Eq.(19:) shows that the mean value of the mode cqi is 



zero. It is the consequence of the fact that the ground 
state of the translationally invariant system is also invari- 
ant under discrete translation |20j . On the other hand, 
the linearized QLEs for the fluctuating operators take the 
followinf forms 



5a 
5c w 

Scqi 



(iAd — n)5a — —G(5c 



+ <5c 10 )+£, (20a) 



/lOj 

(iQoi + 7)^01 - 



-G(5a + 5a^) 



- A t 



01) 



where 



/2 

G = UoVNatf! + -V) 



(20b) 
(20c) 

(21) 



is the effective atom-photon coupling parameter. By 
defining the quadrature fluctuations 

5X a = ^=(5a + 5a^),5P a = ^{5a- 6<J), (22a) 
V2 V2i 

6Xj = -^(5 Cj + 5c]), SPj = -j^Sci + 5c]), (22b) 

with j = (10,01), the QLEs can be written in the com- 
pact matrix form 



= Au(t) +n(t), 



(23) 



where u = [5X a , 5P a , 5Xi , 5P\o, 5X i,5Poi] T is the vec- 
tor of continuous variable fluctuation operators and 



n(t) = fcW,e P W,/.ioW,/pio(t),/xoi(<),/ P oiW] T - 



is the corresponding vector of noises. The 6x6 matrix 
A is the drift matrix given by 
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"01 


~7 
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(25) 



As is seen from the drift matrix, the side mode 
(5X i,5Pqi) has been decoupled from the optical mode. 
It is due to the fact that in the QLEs, Eqs.(18:), there 
is no standard radiation pressure coupling between the 
side mode coi and the optical field. Instead, it is coupled 
to the radiation field via the nonlinear term, — ^-a'acoi 
which is as small as l/y/~N of the radiation pressure term, 
and is deleted during the linearization process (because 
the mean field value of this mode, f3 , is zero). In this way 
the side mode cqi acts as a medium level in the dynamics. 



B. Stationary Quantum Fuctuations 

In order to study the stationary properties of the sys- 
tem it is enough to focus our attention on the subspace 
spanned by the optical mode and the side mode Cio- 
It means that we can consider only the upper block 
of the drift matrix of Eq. (251. The system is stable 



only if the real part of all the eigenvalues of the ma- 
trix A are negative. These stability conditions can be 
obtained by using the Routh-Hurwitz criterion [3U]. Due 
to the linearized dynamics of the fluctuations and since 
all noises are Gaussian the steady state is a zero-mean 
Gaussian state which is fully characterized by the 4x4 
stationary correlation matrix (CM) V , with components 
Vij — (ui(oo)uj(oo) +itj(oo)iti(oo))/2. When the system 
is stable such a CM is given bv|32) 



E 

k,l 



ds / ds'M lk {s)M 3 i{s')D kl { s ~s'), (26) 
o Jo 



where M(s) = exp(As) and Df~i{s — s') — (nfc(s)nj(s') + 
ni(s')rik(s)} /2 is the matrix of stationary noise cor- 
relation functions. For the noise diffusion matrix 
we have Dki(s — s') = Dki5(s — s'), where Dki = 
Diag[/-c, k, r y(2ng + 1), j(2nB + 1)] is the diffusion matrix, 
with fifl = [exp(fijx) m /fcsT) — 1] _1 as the mean num- 
ber of thermal excitations of the Bogoliubov side mode 
(SXio, 5Pio) whose frequency of oscillations is given by 



"10 "io • 



(24) 



Therefore, Eq.(|26|) becomes 

V = I dsM(s)DM T (s). 



(27) 



(28) 
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FIG. 2. (Color online) (a) The mean cavity photon number 
and (b) the mean value fraction of atoms in the Bogoliubov 
side mode Cio versus the normalized cavity detuning A c /ojr 
for two values of uj S w = ur (thin line) and u> sw = Wur (thick 
line). The parameters are N — 6 x 10 4 , Uo = 0.96cj_r, k = 
363.9wh, 7 = 0.001k, r) = 80.06^ and T = lO" 7 ^. 



When the stability conditions are satisfied M(oo) = 
and Eq. ( 28 ) will be equivalent to the following Lyapunov 



equation for the steady-state CM 

AV + VA T = -D. 



(29) 



Equation! 29 ) is linear in V and can be straightforwardly 
solved. However, the explicit form of V is complicate and 
is not reported here. 



IV. NUMERICAL SOLUTIONS TO THE 
MEAN-FIELDS AND FLUCTUATIONS 

In this section we first discuss our results based on the 



numerical solutions of Eqs.(19l for the mean fields and 



then solve Eq. ( 29 ) to obtain fluctuations in the number 
of atoms and photons and their entanglement. We will 
show how the atom-atom interaction affects the mean 
fields, cavity resonance and atom-photon entanglement. 
We analyse our results based on the experimentally fea- 
sible parameters given in Ref. [12 H2] ■ 



A. The mean-field solution 

Here we are going to study the effect of s-wave scat- 
tering on the behaviour of the mean-field values of the 
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FIG. 3. (Color online) (a) The normalized resonance fre- 
quency of the cavity Ai Re!> X<-^R and (b) the normalized effec- 
tive atom- photon coupling G/ujr (thick line) and the mean 
photon number magnified by 10 3 (thin line) versus the nor- 
malized s-wave scattering frequency lj sw /ujr. The parameters 
are the same as those of Fig[2] The cavity detuning has been 
set at A c = 28966cj fl . 



optical field and the Bogoliubov side mode Cio, given 
by Eqs.(19l, quantitatively. We consider N = 6 x 10 4 



atoms distributed in the optical cavity of length L=178 
fim. with bare frequency w c corresponding to a wavelength 
of A=780 nm. The optical mode is coherently driven at 
rate ?7=80.06u;fl with the recoil frequency for rubidium 
atoms lur — 2ttx 3.57 kHz and the one-atom light shift 
C7 =0.96 w R [El [20]. 

In Fig(2] we have plotted the mean number of photons 
and the fraction of condensate atoms occupying the Bo- 
goliubov mode cio versus the normalized atomic detuning 
A c /ljji at a fixed temperature of T = 10~ 7 K, and for 
two values of s-wave scattering frequencies u) sw . As is 
seen the increasing of the s-wave scattering interaction 
shifts the resonance frequency of the cavity to the lower 
values. Besides, it causes the number of atoms in the 
Bogoliubov mode to decrease i.e., leads to the depletion 
of this mode. This is in good accordance with the results 
obtained from numerical solution of the GPE (Fig.l of 
Ref.[33]). 



From Eq.( 19 a) one can obtain the resonance condition 
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A («-) = ^ / °(l + V2A+/3 1 2 ). (30) 

Such a resonance shift has been obtained in Ref. [T3j with- 
out considering atomic collisions. As was mentioned in 
that reference the BEC acts as a Kerr medium that shifts 
the empty-cavity resonance. However, due to the atomic 
collisions the Bogoliubov side mode cio is depleted to 
other modes i.e., the mean value /?j decreases. So a 
BEC with atom-atom interaction has a different reso- 
nance in comparison to a non-interacting one. For clar- 
ity, we have shown the effect of atom-atom interaction 
on the resonance of the cavity in Figj3^a). For a non 
interacting BEC, i.e., uj sw — 0, the resonance occurs at 
Al Rs - ) =28966 ujr while for very large values of U) sw it 
goes to NU a /2 = 28800^. In Figgb) the effective 
atom-photon coupling parameter G (thick line) , given 
by Eq.(21), and the mean photon number of the cavity 
(thin line) have been plotted versus the normalized s- 
wave scattering frequency. This figure shows that both 
the effective atom-photon coupling parameter and the 
mean number of photons decrease with u) aw . All of these 
reductions are the direct consequence of the resonance 
shift exerted by atomic collisions. In the next subsection 
it will be shown that increasing the strength of atomic 
interaction causes the number of fluctuating photons to 
decrease. So the decrease of the total cavity photon num- 
ber is a sign of atomic collisions. These results coincide 
with those obtained by numerical solution of GPE [53] , 



B. Fluctuations and Entanglement 

After calculating the mean-field values we can obtain 
the elements of the drift matrix A and solve for the steady 
solutions of Eq.(23). As explained before, by solving the 



Lyapunov equation (Eq.(29l) we can obtain the corre- 



lation matrix V which gives us the second-order corre- 
lations of the fluctuations. The correlation matrix cor- 



responding to the upper block of A in Eq.(25l can be 
written as 



V 



A C 
C T B 



(31) 



where A and B represent the correlations of the pho- 
tonic and atomic degrees of freedom respectively, and C 
describes the cross correlations. In this way we can cal- 
culate the incoherent excitation number of photons 



3 n ph = (5a' 5a) 



V u + V 22 -1 



(32) 



and the incoherent excitation number of atoms in the 
Bogoliubov side mode 



5n B 



5c\ Sc 10 



V 33 + 144-1 



(33) 
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FIG. 4. (Color online) The incoherent excitation numbers of 
the photons (a), the incoherent excitation numbers of atoms 
in the Bogoliubov side mode cio (b) and the the entanglement 
between the Bogoliubov mode cio and the optical field (c) 
versus normalized cavity detuning A c /ujr for two different 
values of uj sw = ujr (thin line) and ui sw — IOcjr (thick line). 
All parameters are the same as those of Fig[2] 



On the other hand the bipartite entanglement between 
the atomic and photonic degrees of freedom can also be 
calculated by using the logarithmic negativity |35j 



En — max[0, — ln2?7 ], 



(34) 



1/2 



where r)~ = 2" 1 / 2 E(V) - y/T^Vf - 4det ^ is thc 
lowest symplectic eigenvalue of the partial transpose of 
the 4 x 4 CM, V.with £(V) = det.4 + detB - 2detC. 

In Figs|4] (a) and (b) we have plotted, respectively, the 
incoherent excitation number of photons and atoms ver- 
sus the normalized cavity detuning A c /lur for two val- 
ues of uj sw . As is seen, fluctuations are small far from 
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FIG. 5. (color online) The effects of temperature on the inco- 
herent excitation numbers of photons (thin line) and atoms in 
the Bogoliubov side mode (dashed line) as well as the entan- 
glement between photons and atoms (thick line). The cavity 
detuning has been set at A c = 28900cj.r and the s-wave scat- 
tering frequency has been considered to be uj sw = uir. The 
other parameters are the same as those of Fig[2] For clarity, 
the values of Sn p h and En have been multiplied by 100 and 
1000 respectively. 



resonance but they increase rapidly near the resonance. 
On the other hand increasing atom-atom interaction de- 
creases the fluctuations in the number of photons and 
atoms. Since the total number of photons is the mean 
value a 2 that obtained in previous subsection plus the in- 
coherent excitation number of photons, we can conclude 
that the stronger atom-atom interaction, the weaker op- 
tical field of the cavity. In fact the atom-atom interaction 
makes larger shifts in the cavity resonance frequency and 
consequently reduces the cavity field intensity. Hence the 
decrease of the cavity output provides a direct measure 
of the atom-atom interaction within the condensate. 

in Fig|4] (c) the entanglement between the Bogoliubov 
mode cio and the optical field has been plotted versus 
the normalized cavity detuning A c /lor for two values of 
lo sw . Again, the maximum of entanglement occurs at 
resonance. Besides, an increasing in the s-wave scat- 
tering frequency causes the entanglement of the atoms 
and photons to decrease. As it was mentioned above, an 
increase in atom-atom interaction shifts the resonance 
of the cavity which reduces the number of photons of 
the cavity and also causes the depletion of the Bogoli- 
ubov mode. This reduction in the mean values of photons 
and atoms leads to a reduction in the effective atom- 
photon coupling (Fig|3jb)) which causes the entangle- 
ment of atoms and photons to decrease. 

The thermal effects on the incoherent excitation num- 
bers of photons and atoms in the Bogoliubov side mode 
as well as the entanglement between photons and atoms 
have been illustrated in Figj3J As is seen from the fig- 
ure, increasing the temperature causes the fluctuations 
of the number of photons and atoms to increase and on 
the other hand reduces the atom-photon entanglement. 
Furthermore, the thermal effects for the range of temper- 



atures below Q.I^lK are negligible. 



V. SUSCEPTIBILITY AND POWER 
SPECTRUM OF THE BOGOLIUBOV MODE 

Finally, we are going to obtain the power spectrum of 
the Bogoliubov side mode cio and also drive its effective 
frequency and damping rate. We will show that the cou- 
pling between the cavity field and the Bogoliubov mode 
which behaves like a mechanical mirror, leads to the split- 
ting of the normal mode into two modes (NMS). A similar 
theoretical approach has been done in Ref . |36j for opto- 
mechanical Bose-Hubbard Hamiltonian (OMBH) by ex- 
panding the atomic wave field in terms of the Wannier 
functions which is valid only for weak atom-field nonlin- 
earity [37 . Instead, here we use the atomic field expan- 
sion, Eq.(12), in the momentum space. For this purpose, 



we solve the linearized QLEs for the fluctuations in the 
displacement operator 5Xiq as 

5X w (lj) = x (w)F T (w), (35) 

where Fx{oj) is the Fourier transformation of the fluctu- 
ations in the total force acting on the Bogoliubov mode 
and x( w ) is its susceptibility which are given by, respec- 
tively, 

G(k -iui) ^ r \ , G ^d 



F t {uj) 



xM 



A 2 + ( K - 
7 — iuj 



A 2 + ( K 



(-) 

10 



/xio(w) + /pio(w) 



il eff 



u 2 



(36a) 



(36b) 



iujT eff 

where f2 e / f is the effective frequency of the Bogoliubov 
mode, 



O 2 



2 i 2 

7 +w m 



A d G 2 d-\A 2 



- 2 ) 



(A 2 + K 2 - LU 2 ) 2 + 4U! 2 K 2 

and T e ff is its effective damping rate, 



r e// = 2 7 - 



2nA d G 2 n i ~ ) 



10 

2 ) 2 +Auj 2 k 2 ' 



(37) 



(38) 



To calculate the spectrum, we need the correlation func- 
tions of the optical noise sources in the frequency domain 

(kM&V)) = (£pMO')> = kS(u + u/), (39a) 
(&(w)e P (w')) = (&M&(w')>* = iKS(u + u'U39b) 
where we have assumed n p h = 0. For the atomic noise 
sources in the frequency domain we have 

(/xM/x(<«/)) = 2 7 (n B + \)5{u + ut), (40a) 
(/»/>')) = 2 7 («s + \)S(u> + (J), (40b) 



(f x (w)f v {u'))=i 1 8{w + w'), 



(40c) 
(40d) 
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FIG. 6. (Color online) (a) Normalized effective Bogoliubov 
mechanical frequency Q e ff/Lu m , (b) normalized effective Bo- 
goliubov mechanical damping T B ff/u) m and (c) normalized 
power spectrum of the displacement operator of the Bogoli- 
ubov mode versus the normalized frequency ui/uj m for two 
values of u aw — 80ujr (thin line) and u sw — UOur (thick 
line) and for A c = 28700a; r and n — 72.8lur. The other 
parameters are the same as those of Fig(2] 



FIG. 7. (Color online) (a) Normalized effective Bogoliubov 
mechanical frequency fi e ///w m , (b) normalized effective Bo- 
goliubov mechanical damping T e ff/u> m and (c) normalized 
power spectrum of the displacement operator of the Bogoli- 
ubov mode versus the normalized frequency u/uj m for two 
values of uj aw — 80ujr (thin line) and ui aw — 140ur (thick 
line) and for A c = 28700oj_r and k = 24.3^^. The other 
parameters are the same as those of Fig(2j 



where we have omitted the indices of f x and f p for sim- 
plicity. Here, n B = [exp( ~^ c ^ 1 ) — l]" 1 is the num- 
ber of thermal excitations for the Bogoliubov mode and 
u> m is the mechanical Bogoliubov frequency of oscilla- 
tion given by Eq.(p7|. By using the relation S x (uj) — 
i / du J 'e-^+^ t (SX la (u J )5X 10 (u J ') + SX 10 (uj')SX 10 (u J )) 
we can calculate the power spectrum of the displacement 



operator SXiq(u>) as follows 
1 



— \ x (u})\ 2 [4r,(n B + -)-' 



(")2 



(")2 



2 K G 2 (A 2 



(A 2 , + K 2 - LO 2 ) 2 + AiV 2 K 2 



(41) 



The modification of the frequency of the Bogoliubov ex- 
citations of the condensate due to the radiation pressure 



given by Eq. ( 37 1 is equivalent to the optical spring effect 
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in cavity optomechanical systems with movable mirrors 

In Figs. [6ja) and[7](a) we have plotted the normalized 
effective Bogoliubov mechanical frequency fl e f t /ui m ver- 
sus the normalized frequency Lu/uj m for two values of cav- 
ity damping rates k = 72.8wr (Fig[6](a)) and k = 24.3wi{ 
(Fig{7] (a)), and for two different values of ui sw = 80u}r 
(thin line) and lo sw — IAQlor (thick line). As is seen from 
both figures, a higher two-body interaction makes the 
condensate more robust and the Bogoliubov frequency of 
the condensate does not significantly deviate from uj m . 
Besides, a decrease in the damping rate of cavity causes 
the appearance of the peaks in the curves. 

In Figs|6jb) and[7]jb) the normalized effective Bogoli- 
ubov mechanical damping T e f f/u) m has been plotted ver- 
sus the normalized frequency uj/uj m for two values of cav- 
ity damping rates k = 72. 8w^ (Fig{6](b)) and k = 24.30;^ 
(Fig{7] (b)). In the case of Fig|6jbj where the damping 
rate of the cavity is nearly equal to the effective atom- 
photon coupling constant (k ~ G), the effective damping 
rate T e f / is much lower compared to the case of Figj7](b) 
where G > n. It means that the stronger atom-photon 
coupling, the higher atom loss and hence the higher value 
of the effective damping. An experimental observation 
of this phenomenon of light induced back-action heat- 
ing and consequent loss of atoms has been reported in 
Ref.[38j. It was also found that the atom loss rate is 
increased near the resonance. On the other hand, the 
effect of atomic collisions on the effective damping rate 
has been demonstrated in FigsjHJb) andj^b) which show 
that a stronger atom-atom interaction causes the effective 
damping rate to be increased near the resonance points. 
Therefore increasing the rate of atomic collisions can help 
us in cooling of the Bogoliubov mode of the BEC by the 
radiation pressure. 

Finally, we are going to investigate the influence of 
atomic collisions on the phenomenon of NMS in the 
power spectrum of Bogoliubov displacement operator. 
The NMS is associated with the mixing between the fluc- 
tuation of the cavity field around the steady state and the 
fluctuations of the condensate (Bogoliubov mode) around 
the mean field. The origin of the fluctuations of the cav- 
ity field is the beat of the pump photons with the photons 
scattered from the condensate atoms. The phenomenon 
of NMS is observable whenever the energy exchange be- 
tween the two interacting modes takes place on a time 
scale faster than the decoherence of each mode. 

In Figsjefc) and gc) we have shown the normalized 
power spectrum of the displacement operator of the Bo- 
goliubov mode Sx versus the normalized frequency uj/u m 
for two values of cavity damping rate k — 72.8ojr (Figj6j 
(c)) and k — 24.3wfl (Fig[7] (c)), and for two values of 
lu sw = SOljji (thin line) and Zj sw = 140w/j (thick line). In 
the case of Fig|6jc) where the damping rate of the cavity 
is nearly equal to the effective coupling constant of atoms 
and photons (/t ~ G), there is no splitting in the power 
spectrum (neither for lo sw = 80u>r nor f° r u sw = 140lur) 
while in Fig{7](c) where G > k the NMS is appeared for 



lu sw = 140wr. So in order to observe NMS firstly the 
system should be in the strong coupling regime where 
the atom-photon coupling is larger than the decay rates 
of photons an atoms. When this condition is fulfilled, 
the NMS can be observable with increasing atom-atom 
interaction. Normal mode splitting of a system of large 
number of atoms coupled to a ring cavity has been ob- 
served experimentally in the strong cooperative coupling 
regime [33] • The splitting of the normal mode has been 
observed by increasing the number of atoms which leads 
to the increase in the s-wave scattering frequency. 



VI. CONCLUSIONS 

In conclusion, we have done a theoretical investigation 
on the optomechanical properties of a one-dimensional 
Bose-Einstein condensate inside a driven optical cavity 
considering the effects of atomic collisions. Due to the 
dispersive atom-photon interaction the atoms develop a 
band structure in the optical lattice of the cavity. In 
the limit of weak photon-atom coupling the lowest bands 
n = ±1 can be excited by fluctuations due to the atom- 
light interaction. On the other hand, the atomic col- 
lisions scatter the atoms to the higher band and also 
populate states with non zero quasimomentum. It has 
been shown that there is a nonlinear optomechanical cou- 
pling between the nonzero quasimomentum states with 
the optical field proportional to 1/yN which can be 
disregarded in the Bogoliubov approximation where the 
number of atoms in the BEC mode is very large. There- 
fore, there is no radiation pressure coupling between the 
atomic modes with nonzero quasimomentom and the cav- 
ity field in the Bogoliubov approximation. 

In this way we have obtained a simplified optomechan- 
ical model where the Bogoliubov side mode C10 is coupled 
to the optical field through the radiation pressure term 
and the atom-atom interaction Hamiltonian behaves like 
an OPA for the Bogoliubov mode with the non-linear gain 
parameter u> sw . It has been shown that in the strong cou- 
pling regime where the effective atom-photon coupling is 
greater than the cavity damping rate, an increase in the 
s-wave scattering frequency would lead to the NMS in 
the power spectrum of the Bogoliubov mode. 

Furthermore, the atom-atom interaction causes the de- 
pletion of the Bogoliubov mode and also shifts the cav- 
ity resonance frequency which leads to a decrease in the 
mean number of cavity photons. Besides, it decreases 
fluctuations in the number of photons and atoms. Hence 
the decrease in the cavity output provides a direct mea- 
sure of the atom-atom interaction within the condensate. 
On the other hand, a stronger atomic collision rate causes 
a decrease in the effective atom-photon coupling which 
leads to a decrease in the entanglement between the Bo- 
goliubv mode and the optical field. 
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